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We study a two-dimensional exactly solvable non-Hermitian PT— non- 
symmetric quantum model with real spectrum, which is not amenable to 
separation of variables, by supersymmetrical methods. Here we focus at- 
tention on the property of pseudo-Hermiticity, biorthogonal expansion and 
pseudo-metric operator. To our knowledge this is the first time that pseudo- 
Hermiticity is realized explicitly for a nontrivial two-dimensional case. It is 
shown that the Hamiltonian of the model is not diagonalizable. 
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1 Introduction. 

Supersymmetrical techniques has been successfully applied to two-dimensional Quantum 
Mechanics (see [1] and the review paper [2]). For the Hermitian case several real two- 
dimensional models - Morse potential [3], Poschl- Teller potential [4] and some others [5] - 
were studied by means of two different SUSY methods: SUSY-separation of variables and 
shape invariance. Some partial solutions of the spectral problems were obtained by this 
approach. 
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During the last decade intensive study of Schrodinger equation with complex potentials, 
but with real spectrum, was performed by different methods. The pioneer papers [6] initiated 
investigation of PT— symmetric systems (see also the review papers [7]), and afterwards more 
general class of pseudo-Hermitian models was considered [8]. 

One key tool for the complexification of two-dimensional models with real spectra is 
given by the intertwining relations between partner Hamiltonians with supercharges of second 
order in derivatives. A particular class of models - complex singuleir two-dimensional Morse 
- has been found [9] to satisfy SUSY-pseudo-Hermiticity, i.e. 



where the complex supercharges intervene in the HSUSY deformation of the standard algebra 
of SUSY QM. Only partial knowledge of the spectrum and wave functions for this model 
was obtained. 

More recently another complexification of the real singular Morse model was considered 
[10], which from now on we will call reguleirized complex Morse system. This Hamiltonian 
being also involved in second order SUSY intertwining does not satisfy (1), but we point 
out that it fulfills standard pseudo-Hermiticity as defined in [8]: 



with some invertible positive-definite operator rj. In particular, since in this model the com- 
plexification arises from a complex coordinate shift, which also provides a regularization, 
pseudo-Hermiticity is rather straightforward [11]. Due to the fact that this model turns 
out to be solvable, we can now discuss explicitly the biorthogonal expansion based on the 
eigenfunctions and their complex conjugated. 

The structure of the paper is the following. The main results concerning the partially solv- 
able model with real two-dimensional Morse potential are reproduced in Section 2. In Section 
3 the exactly solvable two-dimensional regularized complex Morse model is analysed. This 
model is not PT— symmetric, but its energy eigenvalues are real. After the description of its 



H^Q+ = Q+H, 



(1) 



(2) 
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spectrum and corresponding eigenf unctions, the biorthogonal basis based on these eigenfunc- 
tions and their complex conjugated is studied. The exphcit construction indicates that the 
Hamiltonian is actually not diagonalizable (concerning one-dimensional non-diagonalizable 
Hamiltonians see, for example, [12] and references therein), and biorthogonal basis has to 
be completed in order to provide the resolution of identity. Finally, according to the well 
known prescriptions, the pseudo-metric operator rj and the corresponding positively definite 
(pseudo)inner product [8], [13] space are displayed explicitly. 



2 The partially solvable real two-dimensional Morse 
potential. 

The pseudo-Hermitian model which we want to study in the next Section originates from 
the following Hermitian Hamiltonians: 
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+ Aa^a^, (3) 

+ 4a^a^ (4) 



where the parameters a, A, ct > are arbitrary real numbers, and x± = Xi±X2. 
For the particular choice of parameter a — —1/2, SUSY intertwining relations 



Q+H = HQ- = Q-R 



with the supercharges: 



= 49+9_^2Q;a_^2Q;coth^^9+ + Q;^coth^ 



- A 



d 



dx 



± 



(5) 



(6) 



link the Hermitian Hamiltonian H(x) from (3) to a partner H(x) of (4) which does not 
contain the first term in the r.h.s of (4) and therefore allows for separation of variables [10]. 



The latter Hamiltonian H{x) is straightforwardly solvable with energies expressed in terms 
of two integer positive numbers: 

K,m = en + Cm + a'; ek= -A[l- ^{k + 1/2)]''; k,n,m ^ 0,1, .. . (7) 

V A 



(10) 



All levels of H with n ^ m are two-fold degenerated and the corresponding wave functions 
can be chosen as symmetric and antisymmetric combinations: 

Km = *n,m + *m,n; Km = *n,m - *m,n, (8) 

where the functions ^n,m were defined as: 

^n,m^ Vn{xi)Vm{x2)- (9) 

In turn, rjk are the standard solutions of the one- dimensional Morse problem and can be 
written in terms of confluent hypergeometric functions: 

(^-d"^ + A(^exp{-2ax) - 2exp{-ax)j^r)n{x) = enVnix); 

r/, = exp(-^)(0^"$(-n,2s, + l;0; ^ = ^ expi-ax); (11) 
\fA 

Sn^- n-l/2>0; e„ = -a's^; n = 0, 1, . . . (12) 

a 

The wave functions of H{x) with the same energies (7) are obtained from (8) acting by 
supercharge Q'^ from (6): 

K^^-Q^Km^- (13) 

The operator has singular coefficient functions, and it is antisymmetric under xi X2- 
The two-fold degeneracy of levels (7) of H under n m is not reproduced, in general, in 
the spectrum of H. While the singularities of at x_ = can be compensated by ^n,m 
^n,m(^)> wave functions ^J^^^ may be nonnormalizable. Up to now the hypergeometric 
functions in expressions for the wave functions (9) did not allow to perform a comprehensive 
analysis of the normahzability of all wave functions, i.e. to prove the exact solvability of the 
model (see details and some examples in [10]). 



3 The regularized complex Morse model. 



In order to avoid the singularities at X- — > 0, which hinder the solvabihty of the model, 
it is useful [10] to perform a suitable complex coordinate shift 

x^x + i5] ^=(5,0) (14) 

with 5 small enough (such that a5 E (0, 7r/2)) in order to remove the singularities from 
the real {xi,X2) plane, preserving the normalizability of the functions r)n{xi) from (11) at 
xi — > — oo. After this complex shift the Hamiltonian has obviously a real spectrum, but the 
analysis of normalizability of wave functions is now essentially simplified by the absence of 
singularities in the supercharges. Complexification of both operators is achieved by the 
same shift (14) in definitions (6). Therefore, their mutual Hermitian conjugacy is replaced 
now by 

Q- = HQ^n- (15) 

— * 

The spectrum of the complexified Hamiltonian H{x+iS), which is still amenable to separation 
of variables, coincides with (7), and all eigenfunctions ^n,m are obtained from (9) by the 
same imaginary shift of x. 

Similarly to the Hermitian case, the intertwining relations (5) lead to the eigenfunctions 
^n,m{^ + iS) of the non-separable non-Hermitian Hamiltonian H{x + id) : 

^ngn^+ ^5) = g+(f + .5>f;^)(f + ^5), (16) 

but now, due to the absence of singularity of Q'^ at x_ — > 0, these wave functions are 
normalizable. The corresponding eigenvalues (see (7)) are two-fold degenerate: one can 
choose symmetric or antisymmetric combinations of ^ri,m- The only exclusions arc the levels 
En,n±i, which arc not degenerate, because antisymmetric functions '^nn±i^ being [10] the 
linear combinations of zero modes of Q~^, are annihilated by Q~^. 

It is known [1] that both Hamiltonians H and H obey the dynamical symmetry properties. 
The fourth order operators it! = Q~Q^ and R — Q'^Q~ commute with H and H, respectively. 
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while they do not mix the degenerate wave functions. For the case of complex potentials 
these operators are not Hermitian because of the relation (15). 

In next Section we will need the eigenvalues r„ m of R for eigenfunctions \E'f^^)(x + 
iS). They can be calculated in terms of "one-dimensional energies" e„, e^^ of (7). Indeed, 
separation of variables in operator H with a = —1/2 gives: 

H{x + i5) = hi{xi + iS) + h2{x2) +a'^; (17) 
hi = -d^ - /i = -dl + /l(^e-2°(^i+*^) - 2e-°(^i+^^)) ; (18) 

= -a2' + /2 = -a2' + A(e-''^^^-2e-"^^); (19) 

The explicit form (6) of the supercharges leads to the following expression: 

R = Q-Q+ = (h2-hi + ^C+C_-C+d_-C_d+^(h2-h, + ^C+C_ + C+d_ + C_d+y (20) 

which for a = — I can be transformed by straightforward calculations to: 

R^{hi-h2f + 2a''{hi + h2) + a\ (21) 

It means that the eigenvalues rn,m are expressed as: 

rn,m = (en " ^mf + 2a^(e„ + tm) (^{m - - ij (^{sm + s^Y - ij , (22) 

where the positive parameters s„ were defined in (12). One can notice that for some integer 
n, m eigenvalues rn,m are not positive (operator R is not Hermitian): r„,„ — a^{l — 4s^) < 
for all n (excluding n = — |]), and rn,n±i — for all values of n. 

In general, besides eigenfunctions of the form (16) some additional normalizable eigen- 
states of H could exist, if they would be annihilated by . or if they would be transformed 
by Q~ into nonnormalizable functions. The second option is excluded due to nonsingular 
form of supercharges. The analysis of zero modes of Q~ is performed analogously to in- 
vestigation in [3] (Subsections 4.3 - 4.5) but up to some appropriate changes in that paper: 
— > Q~] h ^ h eid^. The required set of ^„ - hnear combinations of N zero modes 



"^In particular, it means that one has to use in these calculations a = 1/2. 
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ill] I — 0,1, . . . , N oi Q - is constructed by means of SUSY-separation of variables [3] and 
the similarity transformation with function ^1^2(^2 — '■ 

N 

^n = E^nA, (23) 
fe=0 

where bni are matrix elements of B, which satisfy the matrix equation: 

EB = BC. (24) 
In this equation E is diagonal matrix with elements 

K = c„„ = -2aV(l + s„); n = 0, 1, 2, . . . , iV, (25) 
and C is the triangular matrix [3] with elements Cnk, such that 

N 

H^n = Yl Cnk^k- (26) 
A;=0 

The direct algorithm for calculation of bni in terms of known c„fc was also given in [3] . 

One can notice that the eigenvalues En from (25) for the values n = 1,2, ... ,N coincide 
with the eigenvalues -En-i,n of ^n-i,n = Q^^n-i,n from (16), which were found by using 
intertwining relations. It is necessary now to compare the corresponding eigenfunctions ^„ 
and *„-i,n. 

Since the eigenvalues rn-i,n oi R — Q~Q^ vanish for all n = 1,2,. . .,N, the eigenfunc- 
tions ^n-i,n of H are simultaneously the zero modes of Q~ , and therefore must be linear 
combinations of with some unknown coefficients dnk- 

N 

Q^*n-i,n = E«nA; n^l,2,...,N. (27) 

k=Q 

Acting with the H onto both sides of this relation and subsequently equating coefficients in 
front of 0,1 gives: 

TV 

En-l,nanl = Yl (^nkCkf, U, I = 1, 2, . . . , N . (28) 



fe=0 
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In matrix form this equation coincides with (24) up to replacing bnk by Unk, therefore these 
matrix elements also coincide up to a common constant factor. This analysis shows that 
functions Q+^^f.^^ coincide with ^„ for n = 1,2,... N, and the eigenvalues En-i,n = En 
still are not degenerate. 

There is only one additional eigenstate in the spectrum of H not obtained from inter- 
twining relations. It corresponds to n = 0, i.e. Eq — — 2q;^so(1 + Sq). Its wave function - the 
lowest zero mode of Q~ - reads: 

*o = exp(-^^)(66)^°+'(6 - 6)-^ (29) 

Thus, the spectrum of the Hamiltonian H{x+i5) is known: it includes two-fold degenerate 
levels En^m with m 7^ n ± 1, non-degenerate levels E^-i^n with n = 1,2, ... ,N and one 
additional level with energy £"0. 



4 Biorthogonal basis and pseudo-Hermit icity. 

The wave functions ^^^^\x + iS) of H{x + iS) (with separation of variables) and their 
complex conjugate functions (^^^^-'(x + id))* form the so called biorthogonal basis for the 
non-Hermitian Hamiltonian H. The corresponding biorthogonality relations 

= j dXi(f)n{Xi+iS)(f)n'{Xl + iS) j dX2^n{x2)^n'{x2) ^ Snn'^mm' (30) 

can be checked straightforwardly and by comparing the integral along the line Xi G (— oo -|- 
i5, +00 + iS) with the integral along the real axis (with no singularities of integrand between 
these lines). 

The construction of the bound-states-biorthogonal basis by means of the wave functions 
*n,m(^ + i^) and *o(^ + i^) of H{x + i5) together with (^^n,m(x + i5)J and (*o(^ + iS))* 
is much less simple. 
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Due to the property (15), for the complex model the scalar products analogous to (30) 
can be written as: 

( (*n,m(^ + iS)^ I ^n',m'(x + i5)) = 

— * — * 

= (*n,m(^ + I Q~Q '^n',m'{x + iS)) = rn,mK,n'Sm,m' ■ (31) 

In the last equality we used the fact that wave functions \E'n,m arc the common eigenfunctions 
both of the Hamiltonian H and of its symmetry operator R — Q~Q^ with the real eigenvalues 

For all pairs n, m with m 7^ n ± 1 the functions in (31) can be made orthonormal by 
suitable normalization factors, real or imaginary depending on the sign of r„ „j. In particular, 

for m = n one can choose i\rn,n\~^^'^^ n,n{x+i5) and ^i|rn,„|~^/^^„,„(x+i(5)^ as the elements 
of biorthogonal basis (r„^„ < 0). 

No analogous simple prescription works for the functions ^n,n±i{.^ + i^)- The zero value 

of 

fn,n±i-i i-G- the zero value of the integral I i^n,n±i^-'^ d?x, definitely signals incom- 

pleteness of the resolution of identity in terms of (nondegenerate) vectors ^n,m, ^n,m- 
order to give a physical interpretation to the model, one should complete the biorthogonal 
basis by suitable additional vectors. 

Recently the problem of investigation of non-diagonalizable Hamiltonians in one- 
dimensional Quantum Mechanics with non-Hermitian Hamiltonians was studied in papers 
[12] (see also the monographs [14]). Up to our knowledge, not much is known about two- 
dimensional non-diagonalizable Hamiltonians. One can conjecture that the procedure to 
complete the basis is somehow similar to the one-dimensional case. Then, in addition, 
one should consider (in the simplest case) the so called (first order) associated functions 

— * 

^n-i,n{x + iS), which solve the inhomogeneous equation: 

{H - En,n±l)^n,n±l = *n,n±l, (32) 

where the function in r.h.s. is the normalizable eigenfunction of H with eigenvalue En,n±i- 
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Then, due to the second Green's identity (Ostrogradsky-Gauss theorem), the equahties 
{d/dN - normal derivative) 

2 



= j(^n,n±l{x + i5)^ (fx = J({H- En,n±l)^n,n±l^^n,n±l(fx = 
+ f^{^n,n±l-^Kn±l), (33) 

demonstrate that the integral over the large contour in the r.h.s. must be zero for arbitrary 
solution $n,n±i oi (32), irrespectively of the fact that it is normalizable or not normalizable. 

In one dimensional models with discrete spectrum (see [12]) for the normalizable case 
one can complete the biorthogonal basis with ^n,n±i, ^n,n±i with corresponding non-diagonal 
terms in the resolution of identity. Then the Hamiltonian H includes Jordan blocks, and it 
is called non-diagonalizable. 

In the two-dimensional case with discrete spectrum, the general discussion is rather 
complicated. So, we restrict ourselves to the simplest case n = 0, m = 1 in order to provide 
some analytical insight without ambition to propose general theorems. 

In this particular case: 

{H - Eo,i)*Ji = ^Ji; ^0,1 = -2a'so{so - 1), (34) 
where the Hamiltonian with separation of variables is: 

H = -a' {^fd! + e^dl + 6^1 + 6^2 - + ^2^) + («o + ^)(6 + 6) - l) ; (35) 

i\ = exp \-a\xx + i6)\\ ^2 = exp [-ax2), 

a a 

and the wave function reads: 

= exp[-^i±^](66)^°(^ + ^ - (36) 



It is convenient to look for the particular solution $0,1 in the following form: 

= exp l-^-^m^riM^,) + 02(6)), (37) 
10 



where use has been made of separation of variables in Eq.(34). Correspondingly, one obtains 
that the function 0i (and similarly for 02 ) fuUfils an inhomogeneous ordinary differential 
equation: 

(CiVi' - CiV'i + (250 - l)</>i) = ^ - (38) 

The general solution of this equation can be expressed in terms of two linearly independent 
solutions y(^i) and -^(^i) with Wronskian W: 

«Vi(6) = /^y(6)+-(6)+^(6) C dr '^'%(^'^ -y(^^) r^"^^W(^- 

The analysis of asymptotic behaviour of 0i and 02 leads to the conclusion that the function 
is not normalizable and, in addition, the large contour integral does not vanish. This 
is expected since the integral of (^o,i)^ is different from zero, the biorthogonal basis (30) is 
complete, and the Hamiltonian H is diagonalizable. 

Coming finally to the partner model with the Hamiltonian H, we remind that the integral 
of (^o,i)^ is zero. The partner (formal) associated function $o,i = Q'^^o,i turned out to 
be also nonnormalizable^, however the large contour integral vanishes as required by (33). 
Therefore, the problem of completing of the resolution of identity remains open. 

The states \E'o and (^o)* niust also be included in the biorthogonal basis and the resolution 
of identity. It is easy to show that the state is orthogonal to (^n,m)* '■ 

= {%^Jx + i5) I Q-^o{x + iS)) = 0. (40) 

It is difficult to find an analytic expression for the pseudo-norm ((^o)* I ^o) of the state 
(29) but numerical evaluations performed with positive vahics for the parameters Sq and 6 
varying in some limited range indicate that the pseudo-norm does not vanish. 

Summarizing, we have found that the biorthogonal expansion related to Eq.(31) for H 
is incomplete with appearance of states ^n,n±i of zero pseudo-norm. In one-dimensional 
Quantum Mechanics this is associated to non-diagonalizability. In our two-dimensional case 

^It is necessary to remind here again that $o,i is only a particular solution of (32). 
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we have not prooven the existence of associated functions which are normahzable. Irrespec- 
tively of that we have discovered an additional state for H constructed from zero modes of 
Q~ which is pseudo-orthogonal to the other states. This vector definitely should also enter 
the biorthogonal expansion for H. 

Continuing the discussion of pseudo-Hermiticity, an imaginary coordinate shift generates 
this property for H{x + iS), since the following equation holds: 

H(x + i5) ^ H*{x-i6) = exp{-2i6d)H\x + i6)exp{+2i6d). (41) 

Comparing with (2), one can conclude that the explicit form of operator 77 in (2) can be 
written as: 

r]s = exp (+2i5d) = O^O; O = exp (+i5d) = Ol (42) 

In terms of 77 (from now on the dependence on 6 is not written explicitly), the new 
(pseudo)inner product is defined [8] as: 

{n{x)\T{x))r,^{n{x)\vT{x)). (43) 

The precise form (42) of 77 gives for arbitrary fl{x) and T{x): 

{n{x)\T{x))^ = {n{x)\7]T{x)) = {On{x)\OT{x)) = J d^x(n{x + i6)yT{x + i6). (44) 

It is clear now, why the pseudometric 77 is positively definite: for Q = F the 77— norm is equal 
to the integral of \Q{x + iS)\'^. 

5 Conclusions. 

Higher order (nonlinear) SUSY algebra has allowed us to construct an exactly solvable 
two-dimensional non-Hermitian quantum model. We stress that this model is not amenable 
to separation of variables, and it can be considered as a specific PT— non-symmetric com- 
plexified version of generalized two-dimensional Morse model with additional sinh~^ term. 
The spectrum of the model is real. Here we focused attention on the property of pseudo- 
Hermiticity of the model. To our knowledge this is the first time that pseudo-Hermiticity is 
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realized explicitly for a nontrivial two-dimensional case. Following the general results, we also 
studied the biorthogonal expansion and the metric operator associated to pseudo-Hermiticity. 
In particular, it was shown that the Hamiltonian of the model is not diagonalizable. 
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